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In this tutorial we review the basic building blocks of Quantum Information Processing with cold trapped 
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which form a universal set of quantum gates. Different ion qubit choices and their respective gate implemen- 
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5^ 1. Introduction 

I Quantum information processing (QIP) requires the coherent manipulation of the joint quantum 

state of a register of N qubits, in a 2^ dimensional Hilbert space. Some of the resulting states are 
highly entangled superpositions in which the superposition parts represent states that differ in 
many of the qubits. For useful computations N will be mesoscopically or even macroscopically 
large. A superposition of two states that are macroscopically distinct was never experimentally 
observed and is, according to many people, even hard to envisage. The violent contrast between 
the superposition principle, governing the behavior of microscopic systems, and our daily life 
experience at macroscopic scales was already noted in the early days of quantum mechanics [2] • 
The effort of realizing a macroscopic quantum computer is thus also (and primarily) an effort 
T-H to advance the superposition principle onto a macroscopic scale and hence it enjoys the unique 
^ position of having both practical as well as fundamental importance. 

Among the many systems studied as a possible QIP platform, trapped-ions have many advan- 
tages. Atomic ions can be trapped by electric fields in ultra-high vacuum and then laser-cooled 
to extremely low temperatures. Figure [T] shows a linear crystal of eight, trapped and laser-cooled, 
^^Sr"^ ions. Owing to their high degree of separation from any thermal environment, quantum 
superpositions of trapped-ions states have relatively long coherence times. Extremely accurate 
control has been demonstrated over their collective internal states OHj. Trapped ions can be 
laser-cooled to the ground-state of their confining potential, thus providing excellent control of 
their motion as well. Qubits have been encoded into trapped-ion states [3, 4J. In the last decade 
or so, all the basic building blocks needed for coherent control of ion qubits were developed [5]. 
Several quantum algorithms have been demonstrated using up to eight or so ion qubits and, 
from a quantum computing point of view, the challenge that experimentalists in this field are 
currently facing is to find ways to scale this system up to large numbers. 
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The challenge of experimentally realizing a large-scale quantum computer is hard, and it has 
many different fronts. A large-scale trapping architecture has to be developed tools for ap- 
plying parallel quantum gates in multiple trapping regions are needed [7], and so forth. One 
particular challenge is to be able to drive high-fidelity coherent operations on a collective quan- 
tum state of a large number of ion qubits. Obviously, the larger the quantum register is, the 
larger is the space of possible quantum operations. However, the number of possible operations 
in a quantum code has to be finite, and in particular, independent of the number of qubits in the 
register. This problem is analogous to the problem of efficiently compiling a computer program 
into a machine code that implements physical logical operations on a logical register of bits, in 
a classical computer. In QIP, a finite set of unitary operations, which can be concatenated to 
efficiently approximate any operation on an arbitrarily large qubit register, is called a universal 
quantum gate set. This tutorial reviews the basic building blocks of quantum computing oper- 
ations with trapped-ion qubits, with an emphasis on ways to implement a universal quantum 
gate set. 




Figure 1. A linear crystal of eight trapped, and laser cooled ^^Sr"*" ions. 



1.1. 



Universal gate set 



The perfect coherent control of a quantum state in a 2"^ dimension Hilbert space requires the 
ability, up to a global phase, to unitarily connect any two state vectors in this space. This seems 
to be a horrendous task, considering that any such unitary operator is generated by A^^ — 1 
independent generators. Fortunately, it was shown that any such operation can be approximated 
to an arbitrary precision by concatenating a finite number of operators chosen from a small (i.e. 
of size which is independent of N) set ^ ^ . This small and finite set of operators is called a 
universal set of quantum gates. 

Up to a global phase, a pure quantum state of a single qubit can be visualized as a vector 
pointing at the surface of a sphere called the Bloch sphere, illustrated in Figj2| Any unitary 
single qubit operation is therefore equivalent to a rotation R = R{(3, 4>, 6) of the Bloch vector. 
This rotation is specified by the three Euler angles /3, (j), and 9. The angles —tt/2 < (3 < it/2, 
and < (/) < 27r determine the direction of n = (cos /3 cos i?!>, cos /3 sin</>, sin/3), the axis around 
which the Bloch vector is rotated, and 6 is the rotation angle. In the Hilbert space of a single 
qubit, rotations are represented by SU[2] operators, i.e. 2x2 unitary matrices with a unity 
determinant. These rotations are generated by the three Pauli matrices ax, (Jy and dz- The 
rotation R is therefore represented by the operator. 
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where a is the Pauh operators vector. 
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Figure 2. All pure states of a single qubit can be geometrically described as all the points on the surface of a sphere, called 
the "Bloch sphere". All possible single-qubit operations are therefore equivalent to rotations in three dimensions. 



Deutsch was the first to propose the quantum Toffoli gate [8], a rotation i? of a single target 
qubit which is conditioned on the state of two control qubits, as a universal gate set (see Fig. Isk). 




Here the state of the target qubit undergoes a rotation R if and only if the state of both control 
qubits is in the logical Later, Barenco and co-authors proved that single qubit rotations, R, 
together with two-qubit controlled-not (CNOT) gates, shown diagrammatically in Fig. ^p, are a 
universal quantum gate set as well |9j. A CNOT quantum gate performs a 6 = tt rotation of the 
target qubit around the Bloch sphere x axis, i.e. R = ax, only if the control qubit is in state 
Both of these gate set examples are not finite since R is defined through continuous parameters. 
It has been shown, however, that by restricting to a small set of rotation angles, any unitary 
operation in a 2^ dimension Hilbert space can be approximated to accuracy e by concatenating 
a string of length log'^(A'^/e), where c is a constant, approximately equal to 2, of operators from 
the universal gate sets above [U [TU] . 



Figure 3. Quantum circuits of universal gate sets, (a) The Toffoli gate is a controlled rotation. The state of a target qubit 
is rotated only if the state of both control qubits is in the logical | f). (b) Single-qubit rotations and controlled-not (CNOT). 
In a CNOT gate the state of a target qubit is rotated by 180° around the x-axis, if the state of the control qubit is in the 
logical I t>. 



1.2. QIP with trapped ions 




(a) 



(b) 






Qubits that are encoded in the internal states of laser-cooled and trapped ions are promising 
candidates as a quantum computing platform [3]. Trapped-ions are well isolated; thus the cou- 
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pling of environmental noise to their internal electronic states is weak. The long coherence time 
attainable with atomic superpositions has been long recognized and led to the development of 
atomic time standards [llj. Manipulation of atomic superpositions with electro-magnetic radi- 
ation has been extensively studied [T^. Single atomic-qubit rotations were demonstrated long 
before the emergence of quantum information science (e.g. |13j). 

In a seminal paper Cirac and Zoller proposed a method to implement a CNOT gate between ion 
qubits by using their long-range Coulomb interaction [T3]. In the last few years, different schemes 
for implementing quantum gates that can serve as parts of a universal set were theoretically 
suggested |15nl9j and experimentally demonstrated on trapped-ion systems with high fidelity 

[201 m]. 

In this tutorial we review some of the current methods for implementing a universal quantum 
gate set on trapped ion qubits. Here we will focus on the universal gate set suggested in [9], 
i.e., single qubit rotations and two-qubit CNOT gates. In Section [2] we review different ion qubit 
choices, in Section [3] we outline methods of implementing single qubit rotations on the various 
qubit choices. Finally, in Section[4]we describe a method for implementing a two-qubit entangling 
gate. Here we focus on a two-qubit phase gate that is based on spin-dependent light forces. 



2. The ion qubit 

Quantum information is encoded in the internal electronic levels of trapped ions. Ions with a 
single electron in the valence shell have a relatively simple level structure and are therefore well 
suited for this purpose. The relatively manageable electronic level structure also simplifies the 
laser-cooling of such "alkali-like" ions. Once singly ionized, all the elements with stable isotopes 
in group A II earth alkalies, Be^, Mg"*", Ca^, Sr"*", and Ba"*", are good candidate ions with a 
single electron in the zero angular momentum orbital, S, of the n*'* electronic level. All of these 
ions have filled S and P orbitals at the (n — 1) level and no electrons in the D orbital of the 
(n — 1) level, if they are heavy enough to have one. All the radioactively stable, singly ionized, 
group B II transition metals, Zn"*", Cd"*" and Hg^, qualify, with (n — 1) filled S, P, and D orbitals 
and a single electron at the n S level. Finally, the only singly ionized Lanthanide with a single 
electron in the valence shell is Yb"*", with filled (n — 1) S and P and (n — 2) F orbitals and a 
single electron at the n S level. The trapping and laser-cooling of all ions mentioned above were 
demonstrated [22^30] . 



2 

n 




Figure 4. A diagram of the typical level structure of a singly ionized "alkali-like" ions. The ground electronic state is 
the zero angular momentum orbital, S, of the n"* electronic level. An electric-dipole allowed transition to the P orbital 
(typically in the blue-ultra-violet range) allows for laser-cooling, optical pumping and state-selective fluorescence detection. 
In those ions that are heavy enough to have a D orbital in the (n — 1)*'' electronic level, a narrow electric-quadrupole 
transition allows for the realization of an optical qubit. 

Since it is only singly ionized "alkali-like" ions that are currently used for QIP purposes, 
it is worthwhile reviewing common features in their internal level structure. Figure |4] shows a 
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Tabic 1. A list of atomic constants of several of the Hyperfine-ion qubits. Here I is the nuclear spin, 7 is the natural width of 
the -Pi/2 level, /^hf is the S1/2 hyperfinc splitting, ^3/2 ^re the S'1/2 ~^ -^1/2 ^rid S1/2 ~^ -P3/2 transition wavelengths 

respectively. 

schematic diagram of the first few levels of such "alkali-like" ions. The ground electronic state 
of the valence electron is ns "^81/2, where n is the valence shell principal quantum number and 
S is the zero angular momentum orbital. The electron can be optically excited to the P orbital 
via an electric-dipole transition. Owing to spin-orbit coupling, the P orbital splits into two fine- 
structure levels: np '^Pi/2 and np ^^3/2- Tablep^lists the S'1/2 P1/2 and Sij2 P3/2 transition 
wavelengths and the natural spectral width ofthe P1/2 levels for the various ions. 

Earth Alkaline ions heavy enough to have a D orbital in the (n — 1) shell (Ca^, Sr"^, and 
Ba+, as well as Hg+ and Yb'*'), have the nd '^D^^2 and nd ^-^5/2 fine-structure levels between 
the 5 and the P levels. Since the S and the D orbitals both have even parity, the electric-dipole 
transition is forbidden and the two levels are coupled only through their electric-quadrupole 
moment. The lifetime of the D levels (and hence the natural width of the S — )• transition) is 
therefore long (narrow), and typically of the order of one second (one Hz). The Z) — )■ P transition 
is an electric dipole transition and therefore, once excited to the P level, the electron can decay 
both to the states in the S as well as in the D levels with a fixed branching ratio, /. The S ^ P 
transition is typically used for laser-cooling and other QIP manipulations, and in those ions that 
have low-lying D levels, re-pumping lasers, e.g. with a wavelength that is tuned to the D ^ P 
transition, are necessary to prevent the ion from remaining in the long-lived D manifold. 

In odd isotopes or even isotopes with non-zero nuclear spin, / 7^ 0; all the above levels are 
split into two or more hyperfine manifolds. A small magnetic field usually splits these manifolds 
further into different Zeeman sub- levels. 

A qubit is best encoded in a pair of levels that are as resilient as possible to decay and 
decoherence. Since levels in the P manifold have a very short lifetime (in the few nS range), 
qubit levels are typically encoded in two S1/2 states or a 5" and D superposition. Three typical 
choices for an ion qubit are discussed below. 

2.1. Zeeman qubit 

Encoding a qubit in a pair of S1/2 levels has the advantage of a practically infinite spontaneous 
decay time. In ions where 1 = 0, the only states in the 51/2 level are the two Zeeman states 
corresponding to the valence electron spin pointing parallel and anti-parallel to the external 
magnetic field direction [31] . The energy splitting between the qubit levels is Qs^bB — 28 
MHz/mT. Here Qs is the bound electron 5- factor, /U^ is Bohr magneton and B is the magnetic 
field magnitude. An advantage of choosing this qubit is that since the qubit levels are the only 
levels in the S'1/2 manifold, optical pumping out of the qubit manifold does not occur. One 
obvious drawback of the Zeeman qubit is that since the energy difference between the two 
qubit states depends linearly on the magnetic field, magnetic field noise will cause dephasing. 
Furthermore, at low magnetic fields the frequency separation between the two qubit states is 
smaller than the P level natural spectral width and therefore state selective fiuorescence cannot 
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Tabic 2. The tabic lists the 5^/2 ~^ -^5/2 transition wavelengths, the -D5/2 level lifetime and, /, the branching ratio of P 
S / P ^ D spontaneous decay in various optical ion qubits 



be directly applied for qubit state detection. 



2.2. Hyperfine qubit 

In odd isotopes or even isotopes with non-zero nuclear spin, I/O, the S1/2 is split into (4/ + 2) 
different mp states in two hyperfine manifolds: F = I ±1/2. A qubit can be encoded in any pair of 
hyperfine states |32H36j . The manifolds are split at a low magnetic field by A/jj = (/ + l/2)Ahf, 
where Ahf is the ground- state hyperfine constant, which for all the ions discussed here is in 
the few GHz range. Table [2] lists A^f values for various hyperfine ion qubits. A small magnetic 
field removes the degeneracy between different Zeeman states in each of the two manifolds. The 
Zeeman shift at a low magnetic field is given by gp^BBniF, and at all magnetic field levels 
it can be analytically calculated from the Breit-Rabi formula [37J. A clear advantage of the 
hyperfine qubit choice is that pairs of levels can be found such that their energy separation 
does not depend, to first order, on the magnetic field. The most commonly used example is the 
|-F = / + l/2,m/ = 0) and \F = I — 1/2, m/ = 0) states at zero magnetic field, often referred to 
as the "clock transition" states. Working exactly at zero magnetic field is not feasible since the 
degeneracy between different Zeeman states has to be removed to enable spectroscopic addressing 
of various transitions. Consequently, at a low non-zero field, Bq, the magnetic field dependence 
of the clock transition is first order but small. At certain magnetic fields, large enough such that 
the magnetic dipole energy is comparable with hA^f, other transitions can be found that are 
first order insensitive to magnetic field variations |44j . 



2.3. Optical qubit 

Ions that have a low-lying D level have the possibility of encoding a qubit into levels connected 
by the S ^ D optical transition. Since the lifetime of the D level is typically of the order of a 
second, the fundamental limit to the optical qubit coherence time is long relative to a typical 
operation time. However, in contrast with the two previously discussed qubit choices where the 
local oscillator is in the Radio frequency (RF) or microwave range, here the local oscillator is a 
laser. The coherence time of the qubit is therefore also limited by a finite laser linewidth. A laser 
linewidth comparable to the D level natural spectral width (~ 1 Hz) requires laser fractional 
frequency stability of the order of lO"^'^, a non-trivial task. This difficulty renders the optical 
ion qubit choice less common. In ions with / 7^ both the S and the D levels are hyperfine 
splitted. Since the P1/2 level of ions cannot decay to the -D5/2, an optical qubit that is encoded 
on the Si/2 — ^ D5/2 transition enables state-selective fluorescence on the S1/2 P1/2 transition 



(see 

the branching ratio of P — )• S/P — )• D spontaneous decay for the various ions. 



2.4). Table 2.3 lists the S1/2 — ^ D5/2 transition wavelengths, the -D5/2 level lifetime, and /, 



2.4. Initialization and detection 

All the ion qubit choices mentioned above are initialized to a fiducial state using optical pumping 
techniques [3]. Qubit state detection is typically performed using state-selective fiuorescence 
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[31 [38l [39] . Photons are scattered from a laser beam resonant with a transition connecting one 
of the qubit states to a short hved state. The scattered light is collected on a photo-detector. 
Whether or not the photo-detector measures light indicates the state onto which the ion qubit 
superposition collapses. For optical ion qubits, state-selective fluorescence is straightforward 
since the two qubit levels are optically separated. For example, the ion scatters photons from a 
laser beam that is resonant with the 5^/2 ~ P1/2 transition |4L)H42| . The ion cannot decay from 
the Pi/2 to the -D5/2 level and therefore if the qubit superposition collapses on the -D5/2 state, 
then no photons are detected during the detection period. Errors in detection here will be due 
to background photon counts and the finite spontaneous decay probability of the -D5/2 state 
during the detection period. Hyperfine qubit levels are typically separated by much more than 
the P level's natural width (typically a few GHz as compared with a few tens of MHz; see Table 
[2]). State-selective fluorescence is therefore possible on a cycling transition, i.e., a transition in 
which, due to selection rules, the P level decays only to a single S state [391 US]. Errors in this 
case will be mostly due to off-resonance coupling of the detection laser to unwanted states [H]. 
A fairly large magnetic field is needed in order to separate the two states of a Zeeman qubit by 
much more than the P level natural width. Hence, state-selective fluorescence cannot be used 
directly. Detection is typically performed in Zeeman qubits by shelving one of the qubit states 
on a spectrally distant meta-stable state, e.g., one of the -D5/2 states [55Hi7] . Other detection 
schemes that use quantum logic protocols to enhance detection efficiency have been demonstrated 
as well [IHlllSj. 



3. Single qubit gates 

Single-ion qubit gates have been extensively reviewed (e.g. [3llll[5D|). Here we initially provide a 
general description of single-ion qubit gates followed by a discussion of specific implementations 
for the various qubit choices. 

The qubit levels are coupled by electro-magnetic (e.m.) traveling plane waves, i.e., we assume 
that the ion is placed in the far-field of the e.m. radiation. This assumption does not always hold, 
however, and there have been several proposals to implement ion qubit gates using near-field 
e.m. fields [3l [5TH55] . The qubit evolution is determined by the time dependent Hamiltonian, 

H{t) = Ho + Vit). (2) 

Here Hq is the free Hamiltonian of a spin connected to a ID harmonic oscillator (h.o.). 

Ho = ^hujo&z + hujm{a}a + ^), (3) 

where ojq is the qubit levels frequency separation and Um is the h.o. frequency. The time- 
dependent periodic potential V{t) is induced by coupling to the e.m. wave, 

V{t) = hno{a+ + a-)cos{kx-ujt + (j)). (4) 

The coupling constant Qq is also known as the Rabi frequency and a^^~ are the spin rais- 
ing/lowering operators, 

a-=ax + iay = \l){t\ (5) 
a+ = a^- ioy = I t)(i I- 

The wave-vector, frequency, and phase of the e.m. plane wave are k, w, and (/>, respectively. The 
position operator can be written in terms of the h.o. creation and annihilation operators, o and 



2013 15:56 Contemporary Physics ToolBox 



8 Roee Ozeri 

at, 

kx = kxeq + kxoia^ + a) = kxeq + ri{a^ + a), (6) 

The phase owing to the equihbrium position of the ion, kxeq, can be absorbed in (p. Here k is 
the projection of k along the trap direction, xq = ^JfLjlmLOm is the h.o. ground-state width and 
T] = kxQ is the Lamb-Dicke parameter. Moving to the interaction representation and using the 
rotating- wave approximation (RWA), the interaction Hamiltonian is, 

Hint{t) = hno/2a+ exp L(ae-*^"* + a^e''^'"^)] e*^'^-^*) + H.C.. (7) 

The RWA keeps only terms that vary slowly in time and contribute the most to time evolution 
under Hint- Expanding the exponent in Eq. ^ we obtain terms oscillating at multiples of uJm- 
When 6 = SLOm, with integer s, only one of these terms is in resonance and dominantly contributes 
to the ion qubit time evolution. Here the Rabi frequency is modified to, 

nn,n+s = ^n+s,n = ^oUn + s\ e^''('^+'^t) | n) | = J^o^n+s,n , (8) 

indicating the coupling of the two qubit levels and the n and n + s h.o. levels. The coupling is 
modified by, Dn+s,m the Debye- Waller factor. For the ID harmonic oscillator the Debye- Waller 
factor is analytically calculated to be. 



Dn^s,n = e-^y^\\^L\^[{rj^). (9) 




Three different cases are commonly used. The carrier transition occurs when s = and the 
two qubit levels are coupled without changing the h.o. motional state. The Hamiltonian in the 
interaction representation is given by, 

Hcarrier = ^{a+B^^ + a_e-<^) (10) 

The red sideband (RSB) transition is the transition for which 6 = —ujm {s = — !)• Here, when 
the ion qubit spin state is raised, a single h.o. quantum of motion is annihilated. The interaction 
Hamiltonian for the RSB is, 

Hrsb = ^^^(aa+e^<^ + ata_e-*). (11) 

This Hamiltonian is identical to the well-known Jaynes-Cummings Hamiltonian in quantum- 
optics. The third case is that of the blue-sideband (BSB) transition, where 5 = +ujm {s = +1) 
and the raising of the qubit spin is accompanied by creating a single h.o. quantum, 

Hbsb = Hint = ^^^(ata+e^* + aa.e"'*). (12) 

Note that both the BSB and the RSB transitions entangle the qubit spin with its motional state. 

When the deviation of the ion from its average position is much smaller than the radiation 
wavelength, ryy^ ((at + a)^) <^ 1, the Debye- Waller factor ^ can be expanded using as a small 
parameter. In this regime, known as the Lamb-Dicke regime, the Rabi frequencies of the three 



2013 15:56 Contemporary Physics ToolBox 



Contemporary Physics 9 

cases discussed above can be approximated by, 

~ f)o [1 - (n + l/2)ry2] (13) 
Qn-i,n ^ ^Qy/nr] (14) 
Qn+i,n ^ ^oVn + Irj. (15) 

The carrier Rabi frequency r2„^„ is corrected by (n + l/2)r/^, a correction that is small in the 
Lamb-Dicke regime. The RSB and BSB Rabi frequencies depend on the ion motion to a greater 
extent and are proportional to r]^/n and r]\/n + 1 respectively. Note that the RSB Rabi frequency 
vanishes for n = since no h.o. quanta can be further extracted. 

The dependence of the Rabi frequencies on r/ and n can be understood by considering the 
effect of photon recoil on the overlap between different h.o. wavefunctions. Looking at the overlap 
between a momentum-displaced h.o. ground-state, a Gaussian, and either a ground-state (carrier 
transition) or the first h.o. excited state (BSB transition), shown in Fig[5] Here the wavefunction 
is plotted vs. momentum in units of 27r^/xo. In these units the h.o. ground-state width is unity 
whereas the momentum displacement is r] <^ 1. The overlap between the ground-state and the 
displaced ground-state is 1 — ?7^/2 since the Gaussians are displaced by rj around their maximal 
overlap point. The overlap between the displaced Gaussian and the first excited state increases 
linearly with rj since the Gaussian peak is displaced along the excited state's linear slope. The 
dependence on n results from the fact that the distance between nodes in the h.o. wavefunction, 
and therefore also the width of each node is roughly hy/njxo- The overlap between n and the 
displaced n + 1 wavefunctions would therefore also be proportional to -y/n. The \/n dependence 
can also be thought of as resulting from bosonic amplification. 

To avoid entanglement of the ion qubit spin with its motion, on-resonance carrier transitions 



are used to drive single qubit gates. Here the interaction Hamiltonian, Eq. (10), is time inde- 



pendent and the qubit time evolution in the interaction representation is given by, 

|*(i)>mt = e-*^'-*/^|M/(0))i„i = e-^^^-"|*(0))i„t ^ R{B, 0, 0)|^(0))i„t. (16) 

Here Q = iln,nt and the Rabi vector direction are given by the unit vector n = ( cos (p, sin 0,0) 
pointing in the Bloch sphere equatorial plane. The rotation matrix i?(0, (p, 6) represents the 
spinor rotation operator. 



«(0,</>,( 



cos I —ie^'f' sin | 
-ie"*"^ sin | cos | 



(17) 



The origin of the coupling, represented by VLq, can vary and depends on the ion qubit choice. 
Next, we describe three different single-qubit gate coupling implementations for the different ion 
qubit types. 



3.1. Magnetic dipole coupling 

Zeeman or hyperfine qubits have a separation between the their levels that is in the RF or 
microwave range. The two qubit states represent different directions or magnitudes of a magnetic 
dipole and can therefore be coupled by inducing magnetic dipole transitions. The time-varying 
potential in the presence of an oscillating magnetic field is, 

V{t) = -fi- Bo cos{kx-ojt + (/)). (18) 

The magnetic moment of the electron is /i = fj^sigs^ + SlL + gil), where S, L, and / are the 
electronic spin, the electronic orbital angular momentum, and nuclear spin operators, respec- 
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Figure 5. Momentum transfer in the Lamb-Dicke regime, (a) the overlap between a harmonic-oscillator (h.o.) ground-state 
and the same ground-state displaced by rj, owing to photon recoil, decreases by r;^/2. Here, we use the units where xq, the 
ground-state width equals unity, (b) The overlap between the h.o. first excited state and the recoil displaced ground-state 
increases linearly with rj. 

tively, and gs, gi and gj are the corresponding gr-factors. The Rabi frequency here is given by 
Qq = {I \fi ■ BqI I). For a Zeeman qubit fi = gs/J'B^- The Rabi frequency will maximize to 
= 27r X 28Bq MHz when Bq (in mT) is perpendicular to z, e.g. Bq = -Bqx. Figure [G] shows the 
probability of finding a Zeeman qubit, in a ^'^Sr"'" ion, oscillating vs. time, due to magnetic-dipole 
coupling. Since for different ion qubit choices and quantization magnetic field magnitudes the 




Time lix$ec) 



Figure 6. Magnetic dipole coupling between the two states of a Zeeman qubit in a single trapped **Sr+ ion. (a) A diagram 
of the Zeeman qubit levels, (b) A Rabi-nutation curve. The probability of finding the ion-qubit in the | i~) state oscillates 
vs. time. 



qubit separation frequency varies between a few MHz and a few GHz, the wavelength of Bq, 
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A = l/|k| will vary between 0.01 and 1000 meters, where the typical xq is between 10 and 100 
nm. The Lamb-Dicke parameter in this case will be negligibly small, resulting in a carrier Rabi 
frequency is almost nearly independent of the ion motion. The disadvantage is that sideband 
transitions cannot be efficiently driven. Another disadvantage is that this way single-ion qubit 
addressing is impossible. The fidelity of magnetically driven single-qubit gates is mainly limited 
by the amplitude and phase noise of the magnetic RF field at the position of the ion. 



3.2. Two-photon Raman coupling 

Two-photon Raman transitions employ a laser as an optical "carrier" of the RF qubit separation 
frequency. In this scheme, photons are coherently transferred between two optical modes with 
frequencies that are separated by uq + sum as illustrated in Fig. ([T]). The laser frequencies are 
off-resonance, by a detuning Aj, with optically allowed transitions from the qubit levels to several 
excited levels, Cj. With the qubit choices mentioned above, Cj are typically the different states 
in the P1/2 and levels. Thus, the electric fields of the two Raman modes can be written as, 

Er/fe = er/fe^r/fe COs(k^/b • X - W^/^t + (j)^/^), (19) 

where red and blue (r/b) correspond to the lower and higher frequency fields, respectively, and 
e is a unit vector in the field polarization direction. To fully solve for the ion evolution, the 
amplitudes associated with the multiple levels involved have to be calculated. However, when 
the detuning is significantly larger than the coupling strength, the amplitude associated with the 
optically excited states, e^, reduces as 1/Aj and oscillates at two time scales, Aj, and a slow time 
scale that adiabatically follows the ground-state amplitudes. It is thus possible to "adiabatically 
eliminate" the amplitudes of the Cj states and write the equations of motion only in terms of the 
ground-state amplitudes and their derivatives. The problem thus reduces to an effective two-level 
problem involving only the qubit levels. The effective Rabi frequency here is, 

ErEf, ^ \d- er\ei){ei\d- ei\ I) 
^0 = ^2^ ■ (20) 

i 

The time-dependent potential would then reduce to the form of Eq. Q , with (f) = (f)f, — (pr and 
k = kfe — k,.. The carrier Hamiltonian in the RWA will be, 

^^0 = H,nt = ^Dn,n{^+e'^ + a^e-'^) + (A^ - A^)c7,. (21) 



The first term on the r.h.s is identical to Eq. (10). The second term on the r.h.s is due to the 



difference in the qubit levels' light shifts induced by the off-resonance Raman beams, 

^i/t = i^2^ A- + 1/^^ -K-, • ^^^^ 



The differential light shift can be tuned to zero by carefully choosing the beam polarization and 
detuning. 

The Lamb-Dicke parameter and therefore the ability to impart momentum and the sensitivity 
of the Rabi frequency to the ions' motion depends on the angle between the Raman beams. In 
order to efficiently drive sideband transitions, the angle between the Raman beams needs to be 
as large as possible and the difference in their wave vector should be such that k has a maximal 
projection along the trap axis. For the purpose of driving single qubit rotations it is desirable to 
avoid changes in the Rabi frequency owing to ion motion and therefore often a co-propagating 
Raman beam configuration is chosen. Raman beams can be focused to a size limited by the 
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Figure 7. A diagram showing Raman couphng. The Raman beams are detuned from transitions to excited states and the 
difference in energy between the two beams has to match the energy separation between the two qubit states. In order to 
excite ion motion, the two Raman beams cannot be co-propagating. A common choice is for the beams to cross at right 
angle with the difference in their wave-vectors paraUel to the trap axis, as shown in the figure. 



(optical) diffraction limit, typically of the order of a wavelength, and therefore single- ion qubit 
addressing is in principle possible. 

Similarly to the magnetic dipole transition case, errors are caused by classical noises. Laser 
intensity noise and beam pointing fluctuations result in noise in the Rabi frequency. Phase and 
frequency noise between the Raman frequency difference and the ion qubit frequency separation 
also reduces the gate's fidelity. However, in this case, in addition to classical noises, the quantum 
nature of light will further contribute to the gate error [56] . For laser beams the effect of quantum 
noise on the interaction with an atom is attributed to spontaneous scattering of photons. Ground- 
state coherence is only affected by Raman spontaneous scattering events in which a population 
is transferred between different ground-state levels [57], or when the scattering amplitude from 
the two qubit states are different [58]. The spontaneous Raman scattering probability during a 
single qubit gate contributes to the gate error. This contribution reduces for larger Aj values. 



3.3. Optical quadrupole coupling 

Optical qubits are manipulated with optical fields. Consider an optical ion qubit encoded in 
superpositions of the |S'i/2)"T') and 1-05/2, "i') states. Since the two qubit states have the same 
parity, it is only electric quadrupole transitions, which require an electric field gradient across 
the ion, which couple the two levels. Given an optical electric field. 



E = eEo cos(k • x - wq* + , (23) 
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the electric quadrupole Rabi frequency is given by, 

^0 = ^{Si/2,m\ie • r)(k • r)\D,/^,m'). (24) 

Figure [s] shows response spectrum of an optical qubit in a single ^'^Sr'*' ion to a narrow line width 
laser scanning across its transition frequency. Red and Blue side-bands are apparent on both 
sides of the carrier transition. 



m=3/2 ■ 



m=5/2 



n D 



5/2 



m=-5/2 



m=-3/2 



m=-l/2 



n S 



1/2 




20- 



Laser Fluency 



Figure 8. Optical electric-quadrupole coupling between the two states of an optical qubit in a single trapped °°Sr+ ion. 
(a) A diagram of the Optical qubit levels, (b) A spectral response of the ion to the laser frequency. On both sides of the 
transition carrier are red and blue sidebands corresponding to emission and absorption of longitudinal (z) or radial (r) 
quanta of motion. 



Unlike two-photon Raman transitions, here the Lamb-Dicke parameter cannot be tuned to 
zero. Errors in optical quadrupole single qubit gates will result from classical noises in the 
intensity and frequency of the driving laser, and in addition, from a thermal component in the 
ions' motion. Another error results from the finite probability of the level to decay during 
the gate. However, since the lifetime of the D levels is typically ~ 1 sec and much larger than 
the gate time, this error is typically small. 



4. Two-qubit gates 

Thus far, we considered only single-qubit rotations. In addition to signle-qubit gates, the standard 
universal gate set includes CNOT entangling operations. CNOT gates rotate the state of a target 
qubit around the x axis by 180°, depending on the logical state of a second, control, qubit. 
Lexicographically ordered, a general two-qubit state is written as the four vector, 



a|tt)+/3|n)+7l it)+S\ U) 



(25) 



The CNOT operation is then represented by the matrix. 



"l o' 
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(26) 
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Since the roles of the target and control qubits are different, implementing this gate on an 
array of trapped-ion qubits requires a physical operation that can distinguish between the two. 
An entangling operation that is symmetric with respect to both ions, appears to be an easier 
operation to implement. Rotating the target qubit by a Hadamard rotation, which up to a global 
phase can be composed from two consecutive rotations H = e*'^/^^(0, 0, 7r)^(0, 7r/2, 7r/2), both 
before and after the gate, we obtain the operation. 



10 
10 
000 1 
00 10 



[ii ® H2] 



10 
010 
10 
0-1 



(27) 



Here, a it phase shift is imprinted on the | ||) state with respect to the three other collective 
spin states. Note that since this gate and the CNOT arc connected via single-qubit rotations, 
they arc equivalent in terms of a universal gate set. A subsequent 7r/2 rotation of both qubits 
around the z axis yields. 





"100 " 




"1000" 


[^1 (7r/2, 0, 7r/2) ® ^2(7r/2, 0, 7r/2)] 


10 

001 
00 0-1 




i 00 
00 z 
000 1 



(28) 



Here, a 7r/2 phase is imprinted on anti-parallel spin states with respect to parallel spin states. 

Both examples above are entangling operations that rely on imprinting a phase on certain 
collective spin states. It is important to emphasize that whether a certain entangling operation 
is a phase gate or not depends on the choice of basis in which it is presented. As an example, 
the last operation, written in the x basis is. 



[^1(0, 7r/2, 7r/2) ® ^2(0, 7r/2, 7r/2)] 



10 
i 
00 i 
000 1 



[i?i(0, 7r/2, -7r/2) ® i?2(0, 7r/2, -7r/2)] 



gi7r/4 



(29) 
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a collective spin flip operation. 



u\ u) 



-177/4 



V2 



u)+i\n)), 

U) = I tt)- 



(30) 
(31) 



In fact, most of the proposed implementations of an entangling gate for trapped-ion qubits 
can be thought of as a phase gate in one basis and a collective spin-flip operation in another. 



4.1. Quantum phases in cyclic motion 

Entangling phase gates change certain two-qubit collective spin states only by phase factor 
multiplication. One of those collective spin states undergoes cyclic evolution during the gate, at 
the end of which, it returns to its initial state with an acquired phase. The phase factors that 
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accompany the cyclic evolution of a quantum state were first studied by Berry for adiabatic state 
evolution [59j and later by Aharonov and Anandan for the general case [60J. Following j60j . let 
us consider a general quantum state undergoing a cyclic evolution with a period r, 

|^(r)) = e^'^l^'(O)). (32) 

The state evolves according to the time-dependent Schrodinger equation, 

H{t)\^{t)) =ih{d/dt)\'^{t)). (33) 

We define the "bare" state, 

|$(t)) = e-*^W|*(t)). (34) 
such that /(t) — /(o) = (j). We can now write a differential equation for /, 

| = -i(M/(t)|F(t)|vI/(t))+i(cI>(t)|||ci>(t)). (35) 



The phase accumulated by l^*) at time r is given by integrating Eq.(35), 

<t> = -l l\^{t)\H{m{t))dt + i l\^t)\^^\m)dt ^4>D+l. (36) 



There are two different contributions to the total phase. The first term on the r.h.s in Eq.(36) is 
4>D, the Dynamic phase, i.e., the time integration of the energy expectation value divided by h. 
The second term on the r.h.s is the geometric part of the phase, 7. This part is independent of the 
exact Hamiltonian under which the state evolves but rather on its evolution path in parameter 
space. Assuming that at every time instance the wavefunction can be mapped onto a parameter 
space a, \^{t)) = |$q), then 7 can be written as an integral along a closed a-contour. 



7 = i [\^t)\j^Mt))dt = i£{^a\^a\^a) ' da. (37) 



To illustrate the geometric origin of 7, consider a wavefunction evolving around the closed 
contour in parameter space, as shown in Fig. ([9]). The infinitesimal phase difference, A7 between 
\^a) and \^a+Aa) is defined as, 

e-A^ = (38) 



Taking the log of both sides of Eq.(38) yields, 



iA7 ~ (^>a|V„|$Q) • da. (39) 



The geometric phase is thus calculated by integrating Eq.(39) around the a-contour yielding 
Eq.(37). By defining the vector potential 

A^=i{(j)a\ya\(l)a), (40) 

we can write an effective "magnetic field" , 



(41) 
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Using Stokes theorem we then express the geometric phase as the effective magnetic field flux 
through the parameter space area, S, 

1 = - j j ndS. (42) 

Here n is the unit vector normal to S. The encirclement of a finite area in parameter space is 
therefore a necessary condition for accumulating a geometrical phase in cyclic motion. 




Figure 9. A wavefunction evolving around a closed contour in parameter space. The two circles represent the wavefunction 
at a and a + Aa. 



An ion qubit is a (pseudo-) spin connected to a harmonic oscillator. A phase gate such as (28) 
can be implemented by acquiring either a dynamic or a geometric phase or both, during a cyclic 
motion either in spin configuration space or in the h.o. phase space. Cirac and Zoller's seminal 
1995 proposal for an entangling gate with trapped-ion qubits [14J implements a phase gate where 
the phase is geometric and the cyclic motion is in spin configuration space. A later proposal by 
the same authors uses a dynamic phase acquired during a cyclic motion in phase space |61j . 
In the following section we will concentrate on a class of gates, initially proposed by S0rensen 
and M0lmer [ISJ and independently by Milburn and co-authors [16], which uses spin-dependent 
forces to acquire a sum of geometric and dynamic phases during a cyclic motion in h.o. phase 
space. 



4-2. Driven harmonic oscillator 

A classical harmonic oscillator evolves periodically in phase space under the influence of a force, 
F{t) = FQCOs{ujt), oscillating at a frequency slightly off-resonance from its harmonic frequency, 
a; = ujm — S. Starting at rest, x{t = 0) = p{t = 0) = 0, the position and momentum of the 
oscillator at subsequent times are given by, 

Fo / TYl 

x{t) = 2" [cos(a;t) - cos{uJmt)] (43) 

P{t) = 9" [^m sm{ujmt) - oj sin(wt)] . (44) 

Initially the driving force increases the oscillation amplitude until, at time t = vr/(5, the force and 
the oscillator are out of phase. From this time onward the force dampens the oscillator until, at 
time t = 27r/(5, the oscillator returns to its initial state. To make an analogy with the quantum 
case, we normalize the h.o. position by 2xq and the momentum by pq = H/xq and move to a 
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frame rotating at the oscillator frequency. Neglecting terms that are proportional to d/um we 
arrive at, 

x'{t) = ^ cos(a;„t) - ^ sin(a;^i) = ^ [1 - cos(<5t)] (45) 

P'it) = ^ sin(c^^i) + ^ cos(c^^t) ^ ^ [sin(,5t)] . (46) 
zxq Po 2nd 

Wc define the complex function a(t) = x'{t) + ip'{t). The h.o. oscillator motion in the complex 
plane is along the curve, 

(47, 

The h.o. therefore periodically evolves around a circle in the complex plane, the area of which is 



We now turn to solving the quantum, driven, harmonic oscillator. The quantum Hamiltonian 
of a driven h.o. is, 

H = nwm{a)a + 1/2) + F{t)x = Hq + V{t). (49) 

Note that here we assume that the force is spatially uniform and therefore it does not vary across 
the h.o. wavefunction. Moving to the interaction representation and using the RWA we write, 

ihf^\^) I = Vi\^) I, with, 

V, = ^(^a^e''' + ae-'''). (50) 

Vj does not commute with itself at different times; however, for a short enough At, we can still 
write, 

\^{t + At))i = e~i^'^^\^{t))i = e(^""*+^"*«) !*(*))/ = D{Aa)\^{t))i, (51) 
where D{Aa) is the displacement operator, 



oo 

D{am = \a)=e->\'^^\n). (52) 

n=0 

When operating on the h.o. ground-state the displacement operator forms coherent states of 
motion. Here the infinitesimal displacement is, 

Aa{t) = -^Foxoe'^'At = [Ax/2xo + ip/po] . (53) 

The real and imaginary part of the displacement can be identified with its position and mo- 
mentum parts, respectively. To add two infinitesimal displacements, we use the Baker-Hausdorff 
formula. 



(54) 
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for operators A and B that commute with their commutator, [A, [A,B]] = [B, [A,B]] = 0. We 
thus get, 



The displacement operator can accordingly be re-written as, 

N 
i=l 



m § a* da 



(55) 



(56) 



with, a = X^i^i ^CKj- The total displacement at time, t, is therefore, 



a{t) = ^ ^a{t)dt = -^l Foxoe^'^'dt' = ^(1 " e''') 



(57) 



The quantum wavefunction and the classical h.o. therefore follow an identical trajectory in phase 
space, illustrated in Fig. 10 As in the classical case, after time Tg = 2Tr/5, the h.o. wavefunction 
returns to its initial state. 



P' 




] /" 

1 V 


y 
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Figure 10. An illustration of the quantum driven harmonic-oscillator wavefunction trajectory through phase-space. The 
oscillator is startiong from the ground-state and is off-resonantly driven with a periodic force with amplitude _Fo. The area 
that the oscillator is circulating in phase space is equal to 5 = tt [ ^2h5' ] ^ ' 



The total phase that is accumulated by the end of the gate is. 



/m^a^da=^(^)^ (58) 



To implement the phase gate (28), the force magnitude and detuning are adjusted such that 



(f) = 7r/2. Comparing the phase to the area encircled in phase space during the gate, Eq. (48), 
we obtain, 

7 = 25. (59) 

A proportionality between the phase accumulated per cycle and the area encircled by the wave- 
function in phase space suggests that the phase has a geometric origin. 
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To make the connection with Eq.(36) we calculate the dynamic and geometric contributions 
to (p. The dynamic phase is simply, 



1 



{a{t)\Vi{t)\a{t))dt 



(60) 



2h 



r^'(«V^* + a-^*)=^(^)2 = 20. 
Jo ™ 



The geometric part of the phase is, 



7 = i ® (a| V|a) • da. 



(61) 



Here phase space is the parameter space (a complex plane) and the wavefunction trajectory is 
parameterized by the complex displacement a. We evaluate the integrand along the wavefunction 
path by. 



(a|v|a) • da = {a\a + Aa) — 1. 



(62) 



Coherent states are an over-complete basis of wavefunctions and are therefore not mutually 
orthogonal. The overlap between two coherent states a and /3 is, 



e-5(H'+l/3|')+"'/3. 



(63) 



We therefore arrive at. 



{a\a + Aa) - 1 ~ -{a*Aa - aAa*) = iIm{a*Aa). 



(64) 



leading to the geometric phase, 



7 = —Im d) a*da 



(65) 



The total phase acquired during the gate is consequently the sum of partially canceling dynamic 
and geometric contributions. Although of different origin, the geometric and dynamic parts of 
the phase are proportional to each other. The total phase, therefore, is also proportional to 
the geometric phase and accordingly, to the area encircled in phase space [62]. The fact that (j) 
depends on the wavefunction path in phase-space rather than on dynamic properties of the gate 
Hamiltonian contributes to the phase gate robustness against noise. 



Following Eq. ( 40 ) we can define a vector potential. 



A = i{a\^\a) = [-Im{a), +Re{a),0] . 



(66) 



The effective magnetic field (as in Eq.(|41[)) is then perpendicular to the complex plane, i.e. 
B = i?^z and constant. 



5. 



dRe{a) dlm{a) ^ 
dRe{a) dlm{a) 



(67) 



By using the Stokes theorem, Eq. (42 ), we again find that the geometric phase equals minus twice 



the phase-space area encircled, as in equations ( 59 ) and ( 65 ) . The motion of the h.o. wavefunction 



in phase-space (in the interaction representation) is therefore analogous to the motion of a 
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charged particle in a uniform magnetic field. The acquired geometric phase is therefore analogous 
to the Aharonov-Bohm phase. 



4-3. Spin- dependent forces 



Acquiring a phase is not sufficient realizing of a two-qubit phase gate such as in Eq. (28). The 
phase acquired and therefore also the wavefunction trajectory through phase-space have to be 
dependent on the collective two-qubit spin state. To this end, spin-dependent forces have to be 
used. 

Since light couples differently to different spin states, it can be used to apply spin-dependent 
forces on trapped-ion qubits. Spin-dependent forces can be applied to Zeeman or hyperfine 



ion qubits by using the differential light shift, Eq. (22). Here, a spatial variation in the light 
shift induces different forces on ions in different collective spin states. Ions in a collective spin 
superposition move along position-momentum paths that are entangled with the different spin 
states. Since the force evolution is cyclic, spin and motion are disentangled at the end of the 
gate. Here, since the different collective spin states are eigenstates of cr^ (SDiT;^, the gate is referred 
to as a (Tz gate. 

It is somewhat less trivial using the light shift, due to coupling to a third level, to apply spin- 
dependent forces to an optical ion qubit. This is because for a given transition wavelength, when 
one of the qubit levels is connected to a third level, only that level will be light shifted. Forces on 
the two ion qubits will only depend on one of the qubit states. Another example of an ion qubit 
on which such gates fail are "magnetic field insensitive" qubits. Here, since the magnetic field 
couples primarily to the electron spin, and to a far lesser degree to the nuclear spin, "magnetic 
field insensitive" qubits have equal amplitudes of different electron spin contributions to the two 
qubit states [63|- Light also couples primarily to the electron spin and therefore the difference 
in light shifts is very small and is mainly due to the small difference in detuning between both 
qubit states. These small differential light shifts are insufficient to efficiently drive an entangling 
gate. Entangling phase gates on the two ion qubits mentioned above are possible, however, when 
differential forces are applied on spin states that lie on the Bloch sphere equatorial plane. In the 



measurement basis, the action of these gates is similar to Eq.(30). Since this gate is a phase gate 
in the basis of eigenstates of it,^ (8> d'^f,, where a^f, = cos cpdx + sin cpdy, i.e., (j) is the angle form the 
positive X-axis, the gate is often referred to as a gate. 



4-4- ""z gates 

In the case of magnetic field-sensitive Zeeman or hyperfine qubits, cr^ gates can be driven with a 



pair of beams similar to the single-qubit Raman gate case, Eq.(20). Here, however, the difference 
in frequency between the beams is far from the qubit transition frequency, ojq, and is close to 
resonance with the trap motional frequency, ojb—ujr = ^^m+^- In the RWA, no terms are resonant 
and therefore no direct coupling between the two qubit levels occurs. Instead, the only effect of 
the light field is to shift the two qubit levels by A-|-/^, given in Eq. 22 

To form a spatially varying potential, the two beams are crossed such that the vector difference 
of their wavevectors, Ak = kb — k^, is parallel to the axial direction of the trap. By choosing beam 
polarizations mutually perpendicular and with a magnetic field direction that is perpendicular 
to both one obtains a local polarization of the light field that varies periodically between right 
circular, linear, left circular, a linear polarization perpendicular to the previous, and back to 
right circular, over a period of In/Ak. This configuration is often referred to as the lin _L lin 
configuration, known from Sisyphus cooling. The light potential of each of the qubit levels can 
hence be written as, 

^in{x,t) = J(A^/t,+ + A^/^,_) + J[A^/t,+ - A^/^,_] cos{Ak ■ X - Acjt + A(/>). (68) 
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Here A^/-f- _|_ and Aj^/-f- _ are the light shifts of the I and t qubit levels owing to the presence of 
a right and left circularly polarized field respectively. The force on each of the qubit levels is 
therefore, 

Fi/tix, t) = = — [A^/^,+ - A^/^,_] cos(Afcx - Aut + Ac/)) (69) 

With two equal-mass ions in the trap, there are two normal modes of motion along the trap 
longitudinal axis. The center of mass (CM) mode, where the two ions oscillate in phase, has an 
angular frequency that equals that of a single ion, Ucm = ^m- The second mode is the "stretch" 
(ST) mode, where the two ions oscillate out of phase, has an angular frequency of ujst = V^^m 
[64j . The position operator of each ion can be written as, 

Xi = Xeq,i + (aL + acni)±^ (aj^ + ast), (70) 

where aim {acm) and (ast) are the creation (annihilation) operators of the CM and ST modes, 
respectively. The root mean square of the ground-states' spatial spread of the two normal modes 
are XcmO and Xsto- The -|- (— ) signs correspond to the ion with a positive (negative) equilibrium 
X position. This sign change is due to the ions' opposite direction of motion in the ST mode. 
When exciting the CM mode the total force is the sum of forces on the two ions, and when 
exciting the ST mode the total force is the difference in forces. The treatment of a driven h.o. 
above assumes a uniform force. In the Lamb-Dicke regime we neglect the variation of the force 



69 over the ions' wavefunction spread, thus neglecting wavefunction squeezing effects during the 



gate [65]. The force is accordingly considered as uniform with a value given by Eq.(69) at the 
ions' equilibrium position. 



4-5. (T(f, gates 

As discussed in section 4.3, it is hard to apply spin-dependent forces to certain ion qubit real- 



izations such as optical-qubits or hyperfine clock-transition qubits [63]. It is possible, however, 
to apply forces that will depend on their superpositions. In particular, it is possible to ap- 
ply forces that will be equal in magnitude but opposite in direction on the two qubit states: 
1 + 0) = -i=(| f) + e*"^! I) and | - (/>) = ;^(| t) + e~*''^| i)). Here, spin-dependent forces are 
applied in a basis that lies in the equatorial plane of the initial qubit (ion energy eigenstates) 



Bloch-sphere. As shown in Eq. 30 in the qubit basis this operation will look like a collective spin 
rotation. This kind of gate was initially proposed by S0rensen and M0lmer [T^] and demonstrated 
both on hyperfine as well as optical qubits |2H [66] . 



5. Summary 

In this tutorial we reviewed several of the building blocks necessary for QIP with trapped-ion 
qubits. We focused on ways by which a universal quantum gate set can be implemented, and in 
particular on quantum gates which are driven by electro-magnetic fields in the far-field region. 
Here, in order to impart momentum to a cold ion, optical fields are needed. Use of optical fields 
requires relatively narrow line- width lasers. Since atomic transitions of ions are typically in the 
violet to ultra violet spectral regions, lasers are needed at this wavelength range. The advent 
of Gallium-Nitride diode lasers has significantly lowered the cost of near-UV and violet laser 
systems. In fact, for some ion species considered as quantum information carriers, such as Ca"*" 
or Sr+, all the necessary wavelengths for QIP operations are available either directly, or with 
frequency- doubled, diode laser systems. There are, however, several difficulties and disadvantages 
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in using light fields as control fields for QIP. One difficulty is in integrating light into a large-scale 
ion-trap architecture [7J. Realistically, all the routing (e.g. with optical fibers) and switching will 
have to be multiplexed into the trap array structure. This is a non-trivial task and several initial 
steps have been taken in this direction [70J. Another disadvantage is that, although treated here 
in a classical fashion, a light field (as well as any other control field) is essentially a quantum 
field, featuring quantum noise in its operation on a qubit. Quantum noise in the operation of 
a certain control field will pose a quantum limit to the operation's fidelity. In the operation of 
light fields on an atom, this error can be thought of as originating from spontaneous scattering 
of photons during the gate operation [56j. Optical light fields suffer from this type of error to a 
much larger extent than do control fields in the radio-frequency or microwave frequency range. It 
is important to note that our review is not at all exhaustive. For example, in the last few years, 
alternative methods have been suggested, with a few initial demonstrations, in which a universal 
quantum gate set is performed on trapped-ion qubits using only microwave fields [5TH55] . Here, 
momentum transfer to trapped-ions is achieved with strong magnetic field gradients. 

As mentioned in the introduction to this tutorial, despite experimentally realizing of all the 
basic building blocks for a trapped-ion quantum computer, there are many different difficul- 
ties that must be overcome until a large-scale device is finally realized. One such difficulty lies 
in the way decoherence, or error, scales as a function of the quantum register size. Without 
active error-correction, the error probability in a quantum register after a given time would 
scale exponentially with the register size, thus negating large-scale quantum computing. Fortu- 
nately, it has been shown that by using certain error-correction protocols, and given that the 
error-probability per operation is below a certain threshold, large scale quantum computing is 
rendered fault-tolerant [67t i68j. This means that under these conditions, the error in a quantum 
computation would be independent of the register size or the number of computational steps. 
The fault-tolerance threshold is not cast in stone and depends on many parameters such as 
the exact noise mechanism, the error-correction protocol used, the register architecture and so 
forth. However, under rather general assumptions, this threshold has been estimated to be in the 
range |69] . In the last few years, a large effort has been directed toward perfecting all 
the basic operations on trapped-ion qubits, e.g., initialization, detection, memory error, and one- 
and two-ion gates, and reducing the error in these operations as much as possible. Currently, 
typical errors in the gates described above are also in the 10-^ — 10-^ range [HI H^l El EU [72], 
with good prospects for further improvements. 

There is a broad consensus, we believe, that a very large gap exists between the QIP capabilities 
demonstrated so far with trapped-ion qubits, or any other qubit technology for that matter, and 
realizing a large-scale quantum computer. On the other hand, this gap is also widely believed to 
be purely technological rather than fundamental. If this is the case, with more effort and time, a 
quantum computer will eventually be realized. Failure to eventually realize a quantum computer 
will be glamorous because it will indicate a need to modify the laws of quantum mechanics when 
describing large-scale systems. 
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